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FOREWORD 


The  series  of  reports  which  is  entitled  "  Technical  Studies  in  Cargo  Handling" 
is  primarily  a  working  paper  reporting  on  the  progress  of  research  or  the  com¬ 
pletion  of  a  portion  of  a  larger  Investigation.  This  study  is  being  published  in  a 
tentative  form  in  order  to  dlBseminate  the  information  as  quickly  as  possible 
among  the  several  groups  who  are  currently  working  on  related  problemt'.  This 
paper  may  be  expanded,  modified,  withdrawn,  or  published  as  a  report  in  the 
series  entitled  "An  Engineering  Analysis  of  Cargo  Handling"  or  some  other  form 
at  a  later  date. 

The  work  described  in  this  report  was  carried  out  under  the  supervision 
and  technical  responsibility  of  Russell  R.  O'  Neill,  and  is  part  of  the  program  in 
Cargo  Handling.  The  research  conducted  under  the  sponsorship  of  the  Office  of 
Naval  Research,  Department  of  the  Navy,  was  performed  in  the  Department  of 
Engineering.  University  of  California,  Los  Angeles.  L.  M.  K.  Boelter  is  the 
Chairman  of  the  Department. 


Submitted  in  partial  fulfillment 
of  Contract  No.  Nonr  233(07). 


ABSTRACT 


This  report  describes  a  Monte-Carlo  approach  to  the  calculation  of  delays 
in  the  multi-stage  shuttle  process  by  means  of  SWAC,  a  high-speed  digital  com¬ 
puter.  Several  codes  were  developed  for  SWAC  to  generate  the  random  time 
elements,  and  to  calculate  the  delays  in  the  2nd  stage  for  3-  ,  4-  ,  5-  ,  and 
6 -stage  shuttle  processes.  It  was  found  that  the  2nd  stage  delays  did  not  seem 
to  be  influenced  by  the  item  number  but  were  affected  slightly  by  the  number  of 
stages,  the  delays  tending  to  increase  with  increasing  number  of  stages. 


TABLE  OF  CONTENTS 


PAGE 


INTRODUCTION  1 

I  RECURRENCE  EQUATIONS  2 

II  DESCRIPTION  OF  COMPUTATION  AND  RESULTS  6 

m  OUTLINE  OF  CODE  23 

BIBLIOGRAPHY  28 

APPENDICES 

A.  Code  for  3-Stage  Process  29 

B.  Code  for  4-Stage  Process  31 

C.  Code  for  S-Stage  Process  32 

D.  Code  for  6-Stage  Process  33 


I 


notation 


d^(k) 


6^(k) 
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t^(k) 


f  j(k) 


th 

j  «  f«T.  T  at  P  to  receive  the  k —  item 

The  delay  to  T^  incurred  waiting  for  T._^  at  t 

(unit  of  commodity) 

^  *1  «  T  at  P  when  delivering  the 

The  delay  to  T^  incurred  wait^g  for  T.^^  at  .  w 

item  (unit  of  commodity) 

Probability  density  function  of  t 

Frequency  of  discrete  delay  value  2J  computed  by  SWAC 
The  cumulative  frequency  distribution  of 
Subscript  denoting  the  i—  stage  (link) 

The  k~  Item  (unit  of  commodity)  transported 

The  number  of  stages  (links)  in  a  shuttle  process 

The  ti!!-  node  -  the  juncture  of  the  i-1—  and  stages  (Unks) 

The  time  required  for  T^  to  convey  the  k—  item  (unit  of  commodity) 

from  P,  to  P. 
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1+1 


/  tn  P  after  having  delivered  the 

The  time  for  T^  to  return  from  P^^^  to  alter  n 


k—  item  (unit  of  commodity)  to  T^^^ 

The  shuttle  (transporting  agent)  in  the  i  stage  (link) 


The  above  notation  is  consistent  with  the  notation 
in  Cargo  Handling  -  I  .  For  convenience  the  corresponding 
introduced  in  the  series.  An  Engineering  Analysis  of  Car^ 


in  Technical  Studies 
terms  which  were 
jIandUng.  are  included 


in  parentheses. 
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INTRODUCTION 


The  basic  recurrence  relations  for  the  delays  found  in  the  general  N~ stage 
shuttle  process  have  been  formulated  by  R.  Bellman  in  [  1  ]  •  This  formulation 

considers  only  one  shuttle  operating  in  each  stage  with  no  storage  at  the  inter¬ 
mediate  nodes.  No  storage  refers  to  the  requirement  that  the  items  be  trans¬ 
ferred  directly  from  one  shuttle  to  the  next  shuttle.  Delays  to  the  shuttles  will 
occur  then  if  the  two  shuttles  do  not  arrive  at  the  rode  positions  simultaneously. 

A  consideration  of  the  delays  is  important  to  an  understanding  of  the  effective¬ 
ness  of  a  shuttle  process,  as  is  shown  by  R.  R.  O’Neill  in  [2]  .  An  effective¬ 
ness  of  1.  00  is  attained  in  a  shuttle  prcdess  if  the  delays  are  always  zero.  If 
the  shuttle  process  is  a  stochastic  process,  (1.  c. ,  if  only  the  frequency  distri¬ 
bution  of  the  element  times  are  known)  then  the  delays  may  or  may  not  be  zero 
and  in  general  will  also  have  a  frequency  distribution. 

The  distribution  of  the  delays  might  be  expected  to  be  dependent  on  the  number 
of  stages  in  the  process  and  also  on  the  item  number.  The  delays  to  the  shuttle 
involved  in  transporting  the  k^item  would  be  different  from  the  delays  to  the 
shuttle  involved  in  transporting  the  k+1^  item. 


-  ■ 

Numbers  in  square  brackets  indicate  references  in  the  Bibliography. 
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I  -  RECURRENCE  EQUATIONS 


The  basic  recurrence  relations  for  the  delays  found  in  the  general  N-stage 
shuttle  process  have  been  formulated  by  R.  Bellman  in  [l]  .  This  process 


is  shown  in  Figure  1. 

P  P  P 
1  2  3 


•  • 


P  P  P 

N-1  N  N+1 


T  T  T 
1  2  3 


T 

T, 


N-1  N 

FIGURE  1 

There  are  N+1  positions  or  nodes,  P^  ,  P,,  ,  ....  N  shuttles, 

T,  ,  T-  ,  ....  T.,  .  Each  shuttle  T,  operates  between  the  nodes  P.  and  P,  ,  . 

1  2  N  i  1  ItX 

The  last  shuttle,  T..  ,  deposits  the  items  at  the  terminal  P„  ,  . 

N  N  +  l 

Since  there  is  no  storage  provision  at  the  intermediate  nodes,  the  items  can 
only  be  transferred  directly  from  one  shuttle  to  the  next.  A  shuttle  may  therefore 
experience  two  types  of  delays  which  are  defined  as  follows; 

d^(k)  =  the  delay  to  T^  incurred  waiting  for  T^  ^^  at 


6j(k) 


P^  to  receive  the  k—  item. 

the  delay  to  T^  incurred  waiting  for  at 

Pi+1  delivering  the  k—  item. 


The  process  time  elements  are  defined  as 


th 


t^(k)  =  the  time  required  for  T^  to  convey  the  k —  item  from 

0 

^i  ^i+1  • 

t'  ^(k)  =  the  time  required  for  T^  to  return  from  to  P^ 

after  having  delivered  the  k—  item  to 
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As  given  in  [ij  ,  the  general  recurrence  relations  for  the  delays  are 

(1)  d^(ktl)  =  max  [t,_ t  /  _^(k)  -  t.(k)  -  t'.(k)  -  6  .(k)  +  d._ ^(kfl).  O] 

(2)  6.(k+l)  =  max  [  t.^j(k)  f  t.^^  (k)  -  t.(ktl)  -  tj  (k)  -  d.(k^«l)  +  6.^j(k),  oj 

Since  initially  there  is  a  stockpile  of  items  at  Pj  ,  will  experience  no  delay 
in  receiving  items. 

(3)  dj(k)  =  0,  k  =  1.  2. _ n 

Also,  Tj^  will  experience  no  delay  in  depositing  the  items  at  the  terminal. 

(4)  6^(k)  =  0  .  k  =  1.  2. _ n 

And,  as  the  first  item  is  transported  through  the  process,  every  shuttle  will 
experience  no  delay  in  delivering  this  item  since  initially  all  the  shuttles  are  in 
position  ready  to  receive  this  first  item. 

(5)  6j(l)  =  0  .  i  =  1,  2,  - N 

The  delay  to  any  shuttle  awaiting  the  first  item  is  then  simply  the  summation  of 

* 

the  previous  transport  times. 

i-1 

(6)  d,(l)  =  Z  t.(l)  ,  i  >  2 

1  j=i  j 

The  delays  may  then  be  calculated  for  any  i  or  k  by  use  of  (1)  and  (2)  if  the 
transport  and  return  times  are  known. 

As  shown  in  [  2  ]  ,  the  over-all  effectiveness  of  a  shuttle  process  is  the 

same  as  the  effectiveness  of  any  stage  of  this  process  if  there  is  no  storage 
provision  at  the  intermediate  nodes.  This  is  true  since  the  average  flow  of 
items  over  a  period  of  time  must  be  the  same  through  every  stage.  The  dis¬ 
tribution  of  delays  is  important  to  an  understanding  of  the  effectiveness  of  a 
shuttle  process  (or  a  single  stage  of  this  process),  since  as  given  in  |^2j 
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effectiveness  is  defined  as  the  ratio  of  the  theoretical  average  no-delay  round 
trip  time  of  a  shuttle  to  the  actual  average  round  trip  time  of  a  shuttle.  This 
report  describes  the  behavior  of  the  delays,  ^2^^^ 

for  3-  ,  4-  .  5<*  and  S-stage  processes.  In  each  of  these  processes,  the  behavior 
was  Investigated  of  the  delays  to  the  2nd  stage  shuttle  transporting  items  k  >  2, 

3,  4,  101,  102,  103  and  104. 

The  recurrence  relations  necessary  to  determine  these  delays  for  any  k  may 
be  obtained  from  (1)  and  (2)  for  N  =  3.  4,  5  and  6.  The  particular  relations  used 
in  the  computation  of  delays  are  listed  below  for  each  process. 


3-Stage  Shuttle  Process 


(7) 

d2(k+l)  - 

max 

(8) 

6  2(^+1)  a 

max 

4-Stage  Shuttle  Process 

(9) 

d2(k+l)  = 

max 

(10) 

= 

max 

(11) 

dg(k+l)  = 

max 

(12) 

6g(k+l)  a 

max 

5-Stage  Shuttle  Process 

[tj(k-»-l)  t^(k)  -  t2(k)  -  t2(k)  - 

+  t'^(k)  -  t^Ck+l)  -  i^ik)  -  d2(k+l)  ,  O] 

[t^(k+l)  +  t^(k)  -  t2(k)  -  t2(k)  -  62(k)  ,  o] 

[tg(k)  t2{k)  -  t2(k+l)  -  t2(k)  -  d2(k+l)  +  6g(k)  ,  o] 

[t2(k+l)  +  t2(k)  -  tjCk)  -  tgvk)  -  6g(k)  +  d2(k+l/  ,  0  ] 
[t^{k)  +  t'^(k)  -  tgCk+l)  -  tg(k)  -  dgCk+l)  .  O] 


(13) 

(14) 

(15) 

(16) 
(17) 


d2(k+l) 
6  2(k+l) 

dg(k+l) 

63(k+l) 

d^(k+l) 


max  [t^(k+l)  +  t^(k)  -  t2(k)  -  t2(k)  -  6  2(k)  ,  o] 
max  [tg(k)  +  tg(k)  -  t2(k+l)  -  t2(k)  -  d2(k+l)  +  t^gd^)  ,  o] 

max  [t2(k+l)  +  t2(k)  -  tg(k)  -  tg(k)  -  6g(k)  +  d2(k+l)  ,  o] 

max  [t^(k)  +  t^(k)  -  tg(k+l)  -  tg(k)  -  dg(k+l)  +  6^(k)  ,  O] 

max  [tg(k+l)  +  tg(k)  -  t^(k)  -  t^(k)  -  6^(k)  +  dg(k+l)  ,  0^ 
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(18)  6^(k+l)  «  max 

6-Stage  Shuttle  Process 


(19) 

d2(k<fl)  *  max 

(20) 

62(k-hl)  »  max 

(21) 

dg(k-»-l)  =  max 

(22) 

6g(k-»-l)  »  max 

(23) 

d^(k+l)  =  max 

(24) 

6^(k+l)  =  max 

(25) 

dg(k+l)  a  max 

(26) 

6g(k+l)  =  max 

[tg(k)  f  t|.(k)  -  t^(k-^l)  -  t^(k)  -  d^(kfl).  0  ] 


[t^(k+l)  i-t|(k)-t2(k)-t2(k)-  O] 

[tg(k>  +  tg(k)  -  t2(k+l)  -  t2(k)  -  d2(ktl)  +  6g(k),  o] 
[t2(k+l)  +  tgCk)  -  tg(k)  -  tg(k)  -  6g(k)  t  dgCk  +  l),  O] 
[t^(k)  t  t^(k)  -  tg(k-^l)  -  tgCk)  -  dgCk  +  l)  f  6^(k),  O] 

[tg(k+l)  +  tg(k)  -  t^(k)  -  t^(k)  -  6^(k)  +  dg(k+l).  o] 

[tg(k)  i-  tg(k)  -  t^(k+l)  -  t^(k)  -  d^(k+l)  +  6g(k).  o] 

[t^(k+l)  +  t^{k)  -  tgCk)  -  tg(k)  -  6g(k)  +  d^(k+l),  O] 

[t6(k)  +  tg(k)  -  tg(k+l)  -  t|.(k)  -  dg(k^.l).  O] 


<9 
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II  -  DESCRIPTION  OF  COMPUTATIONS  AND  RESULTS 


Hypotheilg  and  Aasumptioni 

The  effectiveneBe  of  the  2nd  stage  shuttle  depends  only  on  the  total  delay 
experienced  by  this  shuttle,  namely,  d^Ck)  ^  bjijt)  .  Reference  to  relations 
(7)  through  (26)  shows  that  each  delay  d^Ck-t-l)  lies  on  a  higher  level  of  com¬ 
plexity  than  ^^(k)  ,  and  each  higher  than  d^fk-i'l).  Therefore  the  com¬ 

putations  are  performed  in  such  a  manner  that  the  results  may  suggest  the  dis¬ 
tribution  of  each  delay  separately.  Specifically,  200  values  of  each  delay,  d^OO, 
b  ^(k) ,  and  d2(k)  -f  ^2^^^  *  computed  for  items  k  >  2,  3,  4,  101,  102,  103 
and  104  in  the  3-  ,  4-  ,  5-  and  6-stage  shuttle  process. 

I 

All  the  delays  are  computed  using  element  times,  t^(k}  and  t^(k) ,  which  are 
assiuned  to  be  obtained  from  the  following  continuous  frequency  function. 

(27)  f(t)  -  l/4e’^^^^,  t  >  0 

I 

All  the  values  of  t^(k)  and  t^  (k)  are  assumed  to  be  random  variables,  with  func¬ 
tional  form  (27),  independent  of  each  other  for  all  values  of  i  and  k. 

Adaption  of  SWAC 

SWAC,  the  hi(^  speed  digital  computer  at  the  Numerical  Analysis  Institute, 
UCLA,  is  used  to  compute  the  desired  delays.  The  operation  of  SWAC  is  briefly 
described  below  for  the  4-stage  shuttle  process.  The  operation  of  SWAC  for  the 
other  shuttle  processes  is  entirely  similar. 

I)  To  compute  d2(2),  62(2)  ,  and  d2(2)  +  62(2)  according  to  relations 

(9)  through  (12) ,  SWAC  generates  10  random  time  elements  ,  t^(l)  ,  t2(l) , 


ilT-( 
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tjd),  t^d).  1^(2)  and  t^U).  Tha  procaaa  of  f anaratinf 

random  variablaa  by  SWAC  is  daacrlbad  in  [2]  on  pagaa  27  to  32. 

SWAC  womputaa  »  t^(2)  -t  t^d)  -  t^d)  '  t^d)*  aubtracta  6^(1)  which  ia 
zero  in  thia  caaa.  from  ,  and  calculataa  AJi2)  ■  max  [0^2“  ^2^^^  *  ^  ]  * 

In  tha  aame  way.  ^  ^  ■  1^(1)  +  t^d)  -  t2(2)  -  t^d)  and  r ^  *  t^d)  +  t^d)  -  tj(2) 

-  tjd)  ara  computed,  and  62(2)  ■  max[^^-d2(2)  +  A^d)  .  0  ]  #  Aji2)  ■ 
max  [d2(2)  -  6^{l),  o]  .  and  62(2)  »  max  [Yj-<1j(2).  O]  ara  obtained 

from  63(1 )  ■  0. 

2)  To  compute  ^2(3) .  .  and  d2(3)  -f  ^2^^^  *  another  aat  of  10  random 

time  alamanta  tj(3) .  tj(2) .  t2(2) .  t2(2) .  t2(2)  .t2(2) .  tjO) .  t^(2) .  t^(2)  uid 
t2(3)«  ia  required.  Since  t2(2)  and  t2(2)  have  alreadbr  bean  generated  in  tha 
pravloua  act.  theaa  two  alamanta  in  the  new  aat  muat  ba  replaced  by  tha  previoua 
valuaa  in  tha  old  aat.  d2(3)  and  6  2(3)  are  than  computed  in  the  same  manner, 
ualng  tha  pravioualy  obtained  valuaa  of  52(2)  and  6^(2), 

3)  SWAC  repaata  the  ab<r#a  procaduraa  for  k  ■  4.  and  punchaa  out  theaa 
raaulta  on  an  IBM  card.  To  obtain  200  groupa  of  theaa  resulta.  SWAC  repaata 

a  • 

the  entire  procaaa  200  timaa. 

4)  To  compute  d2(k) .  62(1^)*  and  d2(k) 62(k)  for  k>  101.  102.  103.  104. 
SWAC  atarta  with  k  ■  2  aa  before,  continuaa  calculation  up  to  k  ■  104.  at  which  it 
punchea  out  the  anawara  only  f or  k  ■  101.  102.  103.  104  and  reaeta  itaalf  for 
another  aequence  of  computationa  atarting  from  k  ■  2.  The  200  groupa  of  delay 
tlmea  ara  obtained  by  200  rapetitiona  of  theaa  procaaaaa. 

Input  and  Output  Data 

Tha  exponential  frequency  function  (27)  ia  approximated  for  computational 


ty 


P 


■  w  .*• 
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purposes  by  the  histogram  shown  in  Figure  2.  This  histogram  Involves  nine 
groups  of  times  assuming  time  Is  a  discrete  variable,  the  mid-points  of  each 
Inteirval  all  being  odd  Integers.  The  area  of  the  histogram  is  divided  into  80 
equal  units  of  area  so  that  the  time,  t,  associated  with  each  unit  of  area  can 
be  stored  In  the  80  positions  In  the  high  speed  memory  of  SWAC. 

The  200  values  of  each  delay  are  interpreted  as  being  200  random  samples 
of  the  delay.  The  computed  delays  will  only  assume  even  Integers  since  by 
relations  (7)  through  (26)  four  values  of  time  are  always  combined.  The  cumu¬ 
lative  frequency  distribution  is  defined  for  every  type  of  delay  as 

^2m+l  IW  ^2j 

where  f„.  is  the  frequency  of  the  delay  value  2j. 

as 

The  frequencies,  f„.  ,  are  obtained  from  the  200  computed  values  of  delay.  The 
cumulative  frequency,  ^201+1  ’  interpreted  as  the  frequency  or  probability 

of  the  delay  being  less  than  or  equal  to  2m-t-l  .  These  dlecret3  distributioris 
Jump  in  value  only  at  odd  numbered  values  of  delay  and  are-  therefore  approxi¬ 
mations  to  the  actual  continuous  cumulative  distributions. 

Time  Required  for  Coding  and  Computation 

The  coding  for  the  3  stage  shuttle  process  required  about  16  hours.  A 
coding  process  includes  planning  of  subroutines,  drawing  of  a  simplified  flow 
diagram,  writing  out  the  code  proper  cell  by  cell,  punching  out  the  IBM  cards, 
tabulating  and  checking.  The  coding  for  the  4-  ,  5-  and  6-stage  shuttle  processes 
required  about  1 0  hours  for  each  process.  Before  starting  the  actual  calculation. 


II  19  15  17  19 

t 

HISTOOMAM  f  (t) 

4.09  4.0 

9.95  4.0 


13  9  7  9 


STATISTIC 
MEAN,  ^ 

STANDARD  DEVIATION,  O’ 


FIGURE  2  HISTOGRAM  FOR  f(t) 


•’  . . . '"■a •''i  ■. 
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SWAC  It  used  to  check  the  code  tftintt  eitx  over-looked  errors  end  mltpltnned 
•ubroutinee.  This  typ€  of  checking  required  about  1  to  1-1/2  hours  for  each 
process.  The  actual  computation  to  obtain  the  200  groups  of  data  for  k  ■  2,  3 
and  4  took  about  5  minutes  for  each  process.  For  k  ■  101,  102,  103  and  104, 
the  same  operations  required  about  10  to  IS  minutes  for  each  process. 

Results 

The  cumulative  frequency  distributions  as  calculated  by  (28)  are  shown  in 
Figures  3  to  14.  The  distributions  are  given  for  each  type  of  delay,  d^Ck)  , 

^^(k),  and  d^ik)  f  ^^(k),  for  the  3-  ,  4-  ,  5-  and  6-stage  shuttle  processes.  In 
each  figure  only  a  few  representative  values  of  k  are  shown  since,  as  can  be  seen 
from  the  graphs,  the  results  are  apparently  random  in  nature  for  different  values 
of  k. 

The  delays  62^^^  increase  with  increasing  number  of  stages.  The 

delays  d»<k)  tend  to  remain  unchanged  with  change  in  number  of  stages.  And, 
the  total  delays  d,,(k)  6,,(k)  also  tend  to  increase  with  increasing  number  of 

m  m 

stages. 
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FIGURE  4  4-3TA6E  PROCESS  DELAY -d«(k) 
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FIGURE  8  4- STAGE  PROCESS  DELAY  -  82^'^) 
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FIGURE  12  4-STAOE  PROCESS  DELAY 
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m  -  OUTLINE  OF  CODE 


1.  Simplified  Flow  Diagram  For  The  3-Stage  Shuttle  Process 


Sub- 

Successor 

Routine 

Operation  or  Successor  Criterion 

Yes 

No 

A 

1 

Generate  and  store  tj(k+l),  tj(k),  t^Ck),  t2(k+i). 

B 

t^Ck),  tg(k).  tg(k). 

B 

Prepare  the  answers  ^ 

C 

C 

Set  tally  for  k  ;  is  k  =  1  ? 

E 

D 

D 

Replace  the  present  t2(k)  by  the  previous  t2(k+l). 

E 

E 

Calculate  =  t^(k+l)  +  tj(k)  •  tgdc)  -  t2(k)  and 

F 

d2(k+l). 

F 

Calculate  |3  f  *  t2(k+l)  -  ^2^^^ 

Q 

62(k+l). 

G 

Set  aside  the  present  t2(kfl)  for  the  following  tgCk). 

H 

H 

Is  k  <  s  (number  of  skips)? 

A 

I 

I 

Convert  the  prepared  answers  into  the  decimal  system. 

J 

J 

Set  tally  for  m  ;  is  m  =  4? 

K 

A 

K 

Punch  out  the  answers;  reset  everything  for  the  next 

L 

A 

repetition;  set  tally  for  r;  is  r  =  100? 

L 

Halt. 
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2.  Simplified  Flow  Diagram  For  The  4-Staffe  Shuttle  Proceaa 


Sub- 

Routine 


Operation  or  Succegeor  Criterion 
Generate  tj(k+l),  t^(k),  t2(k),  t^Ck+l),  t^Ck),  tj(k), 
t2(k+l),  tg(k),  t^<k),  t^(k). 

Prepare  the  answers  ^^(k),  and  ^2^^^ 

Set  tally  for  k;  is  k  =  1  ? 

Replace  the  present  t2(k)  and  t2(k)  by  the  previous 
t2(kH)  and  t2(k+l). 

I  I 

Calculate  o,  =  t,  (k+l)  +  t,  (k)  -  t  (k)  -  t*(k)  and 
k.  1  1  2  2 

d2(k+l). 

Calculate  0  ^  6 2(^+1) 

autid  dg(k+l). 

f  » 

Calculate  y  .  =  t.(k)  +  t.(k)  -  t-Ck-f-l)  -  t,(k)  and 
k  4  4  2  S 

Set  aside  the  present  i2(k+l)  and  tg(k+l)  for  the 
following  t2(k)  and  t2(k). 

Is  k  <  s  (number  of  skips)? 

rr 

Convert  the  prepared  answers  into  the  decimal 


Successor 
Ves  No 


system. 

Set  tally  for  m;  is  m  =  4? 

Punch  out  the  answers;  reset  everything  for  the 
next  repetition;  set  tally  for  r;  is  r  =  100? 


Halt. 
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•3.  Simplified  Flow  Olafram  For  The  5-Stage  Shuttle  Proceaa 


Sub- 

Routine 


Operation  or  Successor  Criterion 
Generate  t^dc+l),  tj(k).  t^dc),  t2(k+l),  t^Ck),  tg(k). 

tg(ki-l),  t2(k),  t^(k).  t^(k+l),  t^(k),  tg(k).  tg(k). 

Prepare  the  answers  d2(k),  62(k),  and  d2(k) 

62(k). 

Set  tally  for  k;  is  k  =  1  ? 

Replace  the  present  t2(k),  t^Ck),  and  t^{k)  by  the 

previous  t„(k+l),  t-(k+l)  and  t  .(k+1). 

«  o  4 

Calculate  ^  ^(k)  =  tj(k+l)  +  t^{k)  -  t2(k)  -  tgik)  and 
d2(ktl). 

Calculate  ^2^^^^  =  "  t2(k+l)  -  t2(k), 

62(k+l)  and  d2(k  +  l). 

* 

Calculate  ^^ik)  =  t^(k)  +  t^(k)  -  t^(k+l)  -  t^ik), 

6-(k+l)  and  d.(k+l). 

0  4 

Calculate  0.(k)  =  t_(k)  +  t_(k)  -  t.{k+l)  -  t.(k) 

4  5  5  4  4 

and  6,(k+l). 

4 

Set  aside  the  present  t2(ktl),  tg(k+l)  and  t^(kfl)  for 
the  following  t2(k),  t2(k)  and  t^(k). 

Is  k  <  3  (number  of  skips)? 

Convert  the  prepared  answers  into  the  decimal 
system. 

Set  tally  for  m;  is  m  =  4? 


Successor 
Yes  No 


£  D 


A  K 
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Punch  out  the  answers,  reset  everything  for  the 


N  A 


nex*  repetition;  set  tally  for  r,  is  r  = 


100? 


Halt. 


Sub- 

Routine 


simplified  Flow  For  The  fi-Stage  Shuttle  Procegg 

Successor 

Oppratton  or  Successor  Criterion 
Generate  tj(k+l),  tj(k).  tj(k),  t2(k+lK  t2(k),  tjCk), 
tjik+l),  t3(k),  t^(k).  t^(k+l).  t^(k).  t5(k),  tj(k+l), 

tjCk).  tg(k).  t'g(k).  • 

Prepare  the  answers  d2(k),  62(k),  d2(k)  +  {^(k). 

Set  tally  for  k;  is  k  =  1  ? 

Replace  the  present  t2(k), 

previous  tjCk+l),  t3(ktl),  t^(k+l).  t5(k+l).  ^ 

Calculate  P  ^(k)  =  t^lk+l)  t  tj(k>  -  t^lk)  -  t^Ck)  and 

d2(ktl). 

Calculate  P^Ck)  =  t3(k)  +  ^k)  -  t^Cktl)  -  t^lk), 

5  (k+1)  and  d_(k+l). 

2  ^  I 

Calculate  Pjik)  =  t^(k)  +  t'^(k)  -  t3(ktl)  -  t3(k). 


6  d^(ktl). 

Calculate  0/k)  =  t^lk)  +  t'^Ck)  -  t^(k+l)  -  t^(k). 

6  and  d^(ktl). 

Calculate  Pjlk)  =  t^lk)  .  ^k)  -  t^lk+l)  -  t'^Ck)  and 

6  ^(k+1). 

Set  aside  the  present  t2(k+l),  t^ik+D.  t^(k+l)  and 


Yes 


.(k+1). 


Set  aside  the  present  t2(k+l),  t^' 


1 
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tg(k-»-l)  for  the  following  t2(k),  t^(k)  and  tg(k). 

K 

Is  k  <  8  (number  of  skips)? 

A  L 

L 

Convert  the  prepared  answers  into  the  decimal 

M 

system. 

M 

Set  tally  for  m;  is  m  =  4? 

N  A 

N 

Punch  out  the  answers;  reset  everything  for  the 

O  •  A 

next  repetition;  set  tally  for  r;  is  r  =  100? 

O 

Halt. 
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